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Abstract 

The amplitude of the deeply virtual Compton scattering off nucleon is computed 
to the twist-3 accuracy in the Wandzura-Wilczek (WW) approximation. The result is 
presented in the form which can be easily used for analysis of DVCS observables. 

Introduction 

Deeply virtual Compton scattering (DVCS) |L], is the cleanest hard process which is 
sensitive to the skewed parton distributions (SPD) and has been the subject of extensive 
theoretical investigations for a few years. First experimental data became recently available 
(see e.g. [@, H |U) and much more data are expected from JLAB, DESY, and CERN in 
the near future. The leading order (1/Q ) DVCS amplitude owing to factorisation theorems 
[@> 0! H can De expressed in terms of twist-2 skewed parton distributions. However, as 
the typical experimentally accessible Q 2 are by no means large, studies of the power (higher 
twist) corrections to the DVCS amplitude are very important. The leading power corrections 
to the DVCS amplitude are of the order 1/Q only (hence twist-3) and their effect on the 
DVCS observables can be rather sizeable. Also the size of the twist-3 corrections increases 
with the increase of the momentum transfer squared t because they enter the amplitude 
typically as \f—t/Q. Therefore, taking into account the twist-3 corrections is important for 
the problem of continuation of the twist-2 SPD's to t — 0. 

Recently the DVCS amplitude was computed in Refs. [II], [II], [T2] including the terms 
of the order 0(1/Q). The inclusion of such terms is mandatory to ensure the electromagnetic 
gauge invariance of the DVCS amplitude to the order yf^i/Q. At the order 1/Q the DVCS 
amplitude depends on a set of new skewed parton distributions. Recently it was shown in 
Ref. [ll] that in the so-called Wandzura-Wilczek (WW) approximation these new functions 



can be expressed in terms of twist-2 SPD's. 

The WW approximation corresponds to neglecting the nucleon matrix elements of the 
'genuine twist-3' operators of the type i^Gip, i.e. neglecting the nonperturbative quark- 
gluon correlations in the nucleon. In principle this approximation is not justified by any 
small parameter of the theory. However, recent measurements of the polarized structure 



function gxix) [13] show that the WW approximation works pretty well in this case. Also 



the estimates of the 'genuine twist-3' contributions to polarized structure functions gr{x) and 



Hl{x) in the theory of instanton vacuum jR, |T^] (which is a model of non-perturbative quark- 
gluon correlations) showed that these contributions are parametrically suppressed relative 
to the 'kinematical' part of the twist-3 contributions by the packing fraction of instantons 
in the vacuum. 

Generically the new twist-3 SPD's in the WW approximation posses discontinuities at the 
points x = ±£ [IB] which potentially, may lead to divergencies in the DVCS amplitude 



and hence to the violation of the factorisation at the 1/Q accuracy However, it was shown 
in [ 12 , |TEj ] for the case of DVCS on the pion target that these dangerous divergencies are 
cancelled to the 1/Q accuracy. In Ref. [|1| general arguments were given that the cancellation 
of the divergencies should persist also for any target. In the present paper we demonstrate by 
explicit calculations that indeed the DVCS amplitude on the nucleon to the 1/Q accuracy 
in the WW approximation is free of the divergencies. The results can easily be used for 
the analysis of the DVCS observables. Collection of formulae the reader can find in the 
Appendix. 

DVCS amplitude on the nucleon at the twist-3 accuracy 

We shall work in the reference frame where the average nucleon momenta and the virtual 
photon momentum are collinear and have opposite directions. For convenience we introduce 
light-cone vectors n, n* such that 

n • n — 0, n* ■ n* — 0, n • n* = 1. (1) 

We also define the transverse metric and antisymmetric transverse epsilon tensor Q 

(-<T)± = ~<T + + n»n*», = e^pnV . (2) 

For an arbitrary vector a M we will use the shorthand notation 

a + = a^n^ 1 , a_ = a^n* u . (3) 

To the twist-3 accuracy the particle momenta can be written as: 

P = ^(p + p')=n*, A—p' — p— -2£P + A ± , 

q = -2£P + ^n, q ' = q -A = ^n-A ± , (4) 

where p,p' are momenta of the initial and final nucleon and q, q' are momenta of the initial 
and final photon respectively. 

The DVCS amplitude on the nucleon is related to the matrix element of T-product of 
two electromagnetic currents 

= -i f d 4 xe-'^ q+c '" )x/2 (p'\TJ^{x/2)r{-x/2)\p) . (5) 



*The Levi-Civita tensor e^ a p is defined as the totally antisymmetric tensor with £0123 
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To the twist-3 accuracy it has been computed in Refs. JIIJ, [TT[| . We use here the form of the 
amplitude obtained in pj in framework of the parton model []. This form is equivalent to 



that of Ref. |TTJ, and can be written as: 



rpflU 



l r 1 



dx 



-<n± 

pu A k 



(Pq) 



+ 



(Pq) 

(q + 2£P) 1 

(Pq) 



(Pq) 



ieLn>F p (x,t)C-(x,t) 



P U A\ 

w 



{F k (x, OC + (x, - iei p F p (x, £)C-(x, £)} 



{F, ± (x, OC + (x, + iei p F p (x, OC'(x, £)} 



(6) 



where the leading order coefficient functions are 

1 



± 



x — C, + ie x + I; — is 



and skewed distributions F^x, £) and F^(x,£,) are defined in terms of the nonlocal light-cone 
quark operators [jj 



F»(x,0 



FJx,C) 



30 dX 1 1 

^-e-™\p'\^(-\n) 1 ^(--\n)\p) , 



00 rlX 1 1 

^e-"\p'm-Xn) 1 ^(--Xn)\p) 

-00 Z7T Z Z 



(7) 



We do not write the sum over quark flavours to simplify notations, the flavour dependence 
in Eqs. (0) can be restored easily by the substitution: 



f=u,d,s,... 

The amplitude (||) is electromagnetically gauge invariant: 

q^ = (q-A) u T^ = 0, 



(8) 



(9) 



to the accuracy 1/Q 2 - In order to have 'absolute' transversality of the amplitude we keep in 
the expression (BI) terms of the A 2 /Q 2 order, applying the prescription of [T/1 ITS] (—g^ — > 

—g^ rp^r) ^ or ^ ne twist-3 terms in the amplitude. Formally such terms are beyond our 

accuracy and they do not form a complete set of l/Q 2 contributions, but we prefer to work 
with the transverse DVCS amplitude. 

^We change little bit notations of Ref. (TTJ 

"'"The gauge link between points on the light-cone is not shown but always assumed. 
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It is easy to see that the third term in Eq. (|^) corresponds to the contribution of the 
longitudinal polarisation of the virtual photon. Note that the last term in the expression (]6|), 
corresponding to transverse polarisation of the virtual photon, is proportional to (g + 2£P) = 
q' + Aj_. Hence, contracting with the transverse polarisation vector e u (q') of the final real 
photon we obtain 

e u (q')(q + 2ZPy = e v (q , )At. (10) 

Therefore, such term does not contribute to any observable with the accuracy 0(A/Q). The 
same contribution has been obtained for the DVCS on the pion. For the case of the pion 



target in Ref. [0 it was shown by a direct calculation that the structure (q + 2£P) is the 
truncated form of the vector q' and hence such term has zero projection onto the polarisation 
vector of the final photon. 

In what follows we compute the skewed parton distributions (0) in the so-called Wandzura- 
Wilczek (WW) approximation. This means that we neglect the contributions of the quark- 
gluon operators, which define what can be called 'genuine twist-3' part. We will demonstrate 
that at least in this approximation the amplitude possesses the factorisation property. An- 
other point, assuming that the contribution of the quark-gluon correlations us smaller than 
the WW part, one can estimate the effects of the 1/Q suppressed terms to various DVCS 
observables. 

Calculation of the matrix elements in the WW approximation 

The aim of this section is the calculation of the matrix elements (0) in the WW- approximation. 
To obtain the answer we make use of the operator identities derived in [11], [L2j, |16j| on the 
basis of technique developed in ref. [|T^] : 

^(x)^ip(-x) = \j^da |aP /t (e~^ x ^ - e^ xV) ) + [e~ a ^ + e^ xV) ) j 

- ie^k^Vi f 1 due {2u ~ 1){xV) f* dae~ a(xV) ^-^(ax)^ 5 ^(-ax) + ... , (11) 

JO JO OXk 

where a = 1 — a and ellipses stand for the contributions of either twist-3 quark gluon 
operators or twist-4 operators^] which we neglect in the WW approximation. An analogous 



expression for the operator i^(x) / y fl j5ip(—x) can be obtained from Eq. (11) by the replacement 
^ — > ^75 in RHS of this equation. Keeping in mind this trivial replacement we shall present 
all derivations only for the vector operator, the generalisation to axial vector operators is 
trivial. 

The symbol T> denotes the derivative with respect to the total translation: 

1 , (12) 



d - 

V a \^)(tx)T[tx, -tx]ip(-tx)j = — {ip(tx + y)T[tx + y, -tx + y}ip(-tx + y)} 



y^O 



§For simple derivation see Ref. [Q 

' In Ref. [[l6j more general relation is derived in which the twist-4 terms of the kind x 2 A 2 are summed 

up. 
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with a generic Dirac matrix structure T and [x,y] = Pexp[igJ dt (x — y) M y4 At (tx + (1 — t)y)\. 
Note that in the matrix elements the total derivative can be easily converted to the momen- 
tum transfer: 

(p'\V^ijj{tx)T[tx, —tx]if)(—tx)\p) = i{p' — p)^{p'\%[){tx)Y\tx ) —tx]ip(—tx)\p) (13) 

and therefore in the matrix elements one can associate T> with the momentum transfer A. 
For what follows it is convenient to rewrite the identity ([TT|) in the more symmetric form: [J 

X> _ 

^(xh^(-x) = j^^(x)^(-x) 



2 Jo L J 



V 

(xV) 



ip(ax)$ip(—ax) 



+ i £ da {e-*^ - e a ^}te^ k x t -^d k i>(ax)^(-ax) + . . .(14) 

To obtain this expression one has to perform an integration by parts in the first and last 
terms of Eq. flnp . 

Our aim is to calculate the matrix elements {p'\4'{z) r j^{ r )^)ip{— z)\p) in terms of the twist- 
2 symmetrical matrix elements {p'\ip{x)^{^^)ip{— x)\p) which appear in the RHS of Eq. (PI). 
At the first step, the symmetric matrix elements should be parametrised for arbitrary vector 
x which is not light-like. The parametrisations of the twist- 2 matrix elements can be written 
as follows: 

?'fT Q/3 T Art 

(p'mx)M-x)\p) = W)Ti[(P*M(*)] + (( 2i J p ))T 2 [(PxU(x)], (15) 
where we use shorthand notations for: 

^) = 4{k}' ((r)) = u(p',s')ru(p,s), (16) 

where U(p, S) and U(p', S') denote initial and final nucleon spinors, respectively. The ampli- 
tudes T\p[(Px), £(x)] can be parametrised in terms of the twist-2 skewed parton distribution 
H,E of Ref. @: 

T 1 [(Pz),Z] = J^due^'lHM, 

T a [(Pz),t] = f^duj^E{u,0- (17) 
Similar formulae can be written for the axial operator: 

mWM>{~*)\p) = Wls))Ti[{PxU{x)\ + ((^p-))T 2 [(Px),Z(x)\, (18) 



II for simplicity we use notation d/dx^ = 
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with amplitudes T\ t 2 given in terms of twist-2 SPD's H, E as: 

Ti[{Pz),$ = J^due^Hiu^), 

f 2 [(Pz),Z] = J^due^^E^O- (19) 



Now, we have to insert the parametrisations (|i~7|) and (p~9[ ) into the RHS of ( |i~4"D and calculate 
the matrix elements up to accuracy 0(Vx 2 ~) which corresponds to the 0(1/Q) accuracy in 
the DVCS amplitude. To collect contributions with such accuracy we first keep the vector 
x off light-cone and after performing the differentiation with respect to vector x we assign 
to the vector x the value on the light-cone x — > z, z 2 = 0. 

Let us illustrate the main steps of the calculations on the example of matrix element of 
the vector operator. The first term in the RHS of the ( JTS ) can easily be computed with the 
result: 

{A^Mx)^{-x)\p) = -^-(( 1+ ))(T 1 [Pz,Z]+T 2 [Pz,Z]) + 

+ ^r((jj))T2[Pz,Z} + (20) 

where ellipses stand for the twist-4 and higher corrections. To get the above expression we 
used that 

^/(Pz) = 7 X = 7+ , (21) 

and the Gordon identity: 

«^^» = -(^)«^» + ^», (22) 

which follows from the equations of motion for nucleon Dirac spinors: 

(jrf - M)U(p, S) = 0, U{p', S') (/ — M) = . (23) 

Consider now a more complicated matrix element of the second term in RHS of (|14|). 
Using (|12"D, ( |i"5|) and (]22f) we obtain for this term the following expression: 



2 nd term = - f da \e~ ia{Ax) + e ia{Ax) \ [a, - -^(xd)] 
2 Jo 1 > (Ax) 



2 Jo I ' (Ax) 

x {UWA{P*U^ (24) 

Note, we can put x = z only after differentiation! The latter can be done with the help of 
the following formulae: 



^T,[«(Px),e(x) 



(Pz) 



Ti[a{PzU], (25) 
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- 7^(^)] T *H^U(*)]U 



1 Af / 9 .9 



"^--feFiW^)^]- ( 26 ) 



" (Ax) v ~~ ^*^v— ,^v~,ji*-* (p z ) 2 ^ V"9a s <9£, 
Applying these formulae and using simple identity: 

P»m = (Pz)((r))-(Pz)((~?D), (27) 
we obtain for the second term the following expression: 

2 nd term 



I/ 1 

2 Jo 






Af 


((7 + }> ( 


1 + a— 


+ 2« 


A^ 




1+a 7T 


2£ 





(28) 



Using the definition of the amplitudes T 12 Eqs. ([T7|) this term can be represented in the form 
which is convenient for the transition to the momentum space: 



2 nd term = J du G^u, f) J da [e l2{Px)[a{u 



-0+5] + e i2(Px)[a(u+Z)-tn 



(29) 



where 



du dt; 



A' 



2s c 



«7+» 



u _d_ + d_ 
du <9£. 



(30) 



Manipulations with the nucleon matrix element of the third term of the operator identity ([14]) 
are similar. For this term we obtain: 



3 d term = ie^ k 



1 du G k (u, jf * rf« { e «(*-O[«C«-0+«] - e a(^)[«(«+0-S } , ( 3 1) 



where 



i + u A + ^ 

<9u d£. 



2s c 



((7+75)) 



du dt; 



H(u,0- 



(32) 



Collecting all terms (p20|) , (p29|) and ([31]) and performing the Fourier transformation on the 
light-cone we obtain the expression for (see definition ([/])) in the WW approximation in 
terms of twist-2 SPD's H, E and H, E: 

F™ W (x,0 = ^((^}}E(x,0-^((l + })(H + E)(x,0 + 



+ duG„(u,^)W+(x,u,^) + ie ± ^ du G*(u,£)W-(x,u,Z) , (33) 
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where we introduced a convenient notation for the WW-kernels: 

W ± (x,u,0 = -\8(x>0- -^-0(z<£)- 1 r 

2 1 u — t u— t J 



, 1 f/>/ / _,6{u>x) n/ 6{u < x) 

± e j > - 1 J - o i < -o 1 , / 

2 1 w + 4 w + 4 



The function G»{u,£) and G k (u,£) are defined in f^) and fl3|). 
The analogous calculation for the axial matrix elements gives: 



WW, 



A 



<75 



f((^))E(x,0 



A, 
2£ 



(34) 



+ j\uG^u,0W + (x,u^) + ie ±llk j\uG k (u^)W4x,u^). (35) 

Let us note that the spinor structures in the functions G^{u, £) and G h (u,£,) are not in- 
dependent. They can be related to each other with the help of equations of motion for 
Dirac spinors (p3|) . However, we prefer not to reduce the number of structures to retain the 
symmetric form of the results ( |33| , |35| ) . 

The generalisation of the WW relations for the non-forward case was discussed first in 
Our results 



Ref. 20 



(l33|J35| ) do not agree, however, with those of Ref. p0| . The reason for 
the disagreement is that in Ref. [20] instead of the complete operator identity of the type fllTD 
only its part which survives in the forward limit has been considered. This simplification is 
justified in order to reproduce only the WW relations in the forward case. The WW relations 
obtained in Ref. |TTJ can be reduced to the form (|33|j35|) obtained here. 



Properties of the WW kernels 

In the previous section we introduced the notion of the WW-kernels, see Eq. (|3~iD. In this 
section we consider the general properties of these kernels. Let us introduce the following 
notations for the action of the WW kernels on a function f(u,£): 



W±®f[x,£] 



-i 



du W±(x, «,£)/(«,£) , 



(36) 



with W± given by Eq. (p4|). We shall call 'WW transform" the resulting functions W± <S> 
f[x, £] and we shall call the 'WW transformation" the action of the WW kernels. 
Limiting cases. We consider two limiting cases of the WW transformation: the forward limit 
4" —>■ and the 'meson' limit £ — > 1. In the forward limit we easily obtain: 



limWl ® f\x,£\ 
lim W_ ® f\x,£\ 



6(x > 0) 
0. 



1 du 



u 



f(u,£ = 0)-0(x<0) 



c du 



1 u 



/(«,e=o) 



(37) 



We can see that the action of the W+ in the forward limit reproduces the Wandzura-Wilczek 
relation for the spin structure function g? [21]. The term with 9(x > 0) corresponds to the 
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quark distributions and the term with 9{x < 0) to the antiquark ones. The W- kernel 
disappears in the forward limit, so that this kernel is 'genuine non-forward' object. 

In the limit £ — > 1 the skewed parton distributions have properties of meson distribution 
amplitudes. In this limit the WW transforms have the form (we use notation f(u,£ = 1) = 
tp(u)): 

gm^ */[*,«] = \{f^ ufW) ±[^- u *„)}, 

which corresponds to the WW relations for the meson distribution amplitudes derived in 

We can see that WW transforms of skewed parton distributions interpolate between WW 
relations for parton distributions and the meson distribution amplitude. Also general form 
of the WW kernels ( |3~4"D allows to derive WW relations for distribution amplitudes of a meson 
of arbitrary spin. 

Mellin moments: One can easily derive the Mellin moments of the WW transform. The 
result has the form: 



dxx N W±®f[x,£\ 



N + 1 



du 



u N+l 



N+l 



± 



u 



N+l 



-0 



N+l 



/(«,£)• (3? 



From this simple exercise we can see important property of the WW transformation, namely, 
if the function f(u, £) satisfies the polynomiality condition, i.e.: 



rl 

J du u N f(u, £) = polynomial in £ of the order N + 1 , (39) 

its WW transform also satisfies the polynomiality condition. Let us consider now the special 
cases AT = 0, 1 of the Mellin moments. These particular moments of the twist-3 SPD's (^) 
do not receive contribution from the 'genuine twist-3' quark-gluon operators, therefore the 
WW approximation for them gives exact results. In the forward limit this observation leads 
to Burkhard-Cottingham (N = 0) p4f and Efremov-Leader-Teryaev (N = 1) |25j sum rules 
for polarized structure function gj*. The corresponding generalizations of these sum rules in 



the case of skewed parton distributions have been discussed in Refs. [0, |16| . From general 
expression (|38|) one gets: 



J 1 dx ®f[x,£] = 0, 



1 dx W+®f[x,£] = f X duf{u,Z) 



-i J-i 

£ rl 



dxxW.® f[x, e] = J / du f(u, , (40) 
l 2 J-i 

1 1 rl 

dx x W + ®f[x,£] = - duuf(u,£). 
l 2 J-i 

Below we shall use these results for the consistency checks of our results for skewed parton 
distributions and for derivation of the sum rules for the distributions F M and F M (0). 
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Discontinuities. In Ref. [IB it was demonstrated that the twist-3 skewed parton distributions 
in the WW approximation exhibit discontinuities at the points x = ±£. This feature is 
related to the properties of the WW kernels. Let us compute the discontinuities of a WW 
transform at the points x = ±£: 

WW±®m + S^\-W ± ®f[i-5^\\ = -vp -f(u,0, 

8^0 2 J -l U — q 

I du 

li m [W ± ®f[-Z + 5,t]-W ± ®f[-t-5,{)) = ±-vp __/(«,£). (41) 

<5^0 I J -i U + 4 

Here vp means an integral in the sense of valeur principal. We see that for a very wide 
class of functions /(«,£) the discontinuity of the corresponding WW transforms is nonzero. 
This feature of the WW transformation may lead to the violation of the factorisation for the 
twist-3 DVCS amplitude. 



Properties of the skewed distributions F„ and 



In this section we analyse the general properties of the skewed parton distributions F^ 
and F^ w (|^) in the WW approximation. We consider the limiting cases and derive sum 
rules for these SPD's which are valid beyond WW approximation. 

Forward limit: In the forward limit the 'vector' SPD F^ in the WW approximation ( |33|) 
vanishes, because unpolarized parton densities of twist-3 are absent. The forward limit for 
the 'axial' SPD F^ w , see (p5|), is nontrivial: 



\imF^ w (x } = 2S- 



lim F„ 

A^O M 



WW , 



Aq(x) + f 1 —Aq(u) 

Jx U 

Aq(\x\) + f 1 —Aq(u) 

J\x\ U 



for x > , 
for x < . 



(42) 



It is easy to recognize in these equations the Wandzura-Wilczek relations for gx{x) ||21|| . In 
the derivation ( |4"2"D we used that twist-2 skewed distribution H is reduced in the forward 
limit to the polarized quark (for x > 0) and antiquark (for x < 0) densities Aq(x) and 
Aq(x). 

Mellin moments: The lowest N = 0, 1 Mellin moments of the SPD's F^ and F^ (^j have 
been analysed previously in Ref. [ID|. Since for these moments the contributions from the 
twist-3 quark-gluon operators are absent, we have to reproduce sum rules of Ref. JTIJ from 
the WW approximation for the functions F^ and F^, see Eqs. (|33|),(p5|). This would be a 
powerful check of our results. 

As it follows from the definitions of twist-3 SPD's (0), the N = moments are reduced to 
the nucleon matrix elements of the local currents which are parametrized in terms of electric, 
magnetic, axial and pseudoscalar form factors of the nucleon: 



du 

du Fpfat) 



((7,))^i(A 2 ) + ((^^))F 2 (A 2 ), 
((7,7 5 ))G a (A 2 ) + ^:(( 75 ))Gp(A 2 ) 



(43) 
(44) 
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Using the properties of the WW kernels ( f4"0| ) one can easily see that the expressions for F 



WW 



fl33|) and for F^ w (|35|) satisfy the general QCD sum rules (|43|). This is a sensitive check of 
our approach. 

The general QCD expression for the N — 1 moment of the SPD's F^ and can be 
obtained from the WW approximation because the quark-gluon operators do not contribute 
to this moment. Again with help of ( |40| ) we can obtain the results for the N = 1 moments of 
the transverse components of the functions F^ and Fj^ which enter the DVCS amplitude 
for the longitudinally polarized virtual photon. For the function F fl± we have 



duuF™ w (u,0 = «7^»~ 



y 1 du u {H(u, + E(u, 0} + G A (A 2 



(45) 



We see that the 'transverse' SPD F Mx is sensitive to the combination 

J du u (H(u, + E{u, 0) + G A (A 2 ) 
which is related to the spin structure of the nucleon in the forward limit: 
fim ^ du u {H(u, + E(u, 0} + G A (A 2 ) = 2J q + Aq , 



(46) 



where J q is a fraction of the total angular momentum of the nucleon carried by quarks and 
Aq is the corresponding fraction of the spin. To derive Eq. (^) we made use of Ji's sum 
rule 0. Also the new angular momentum sum rule (|4"5| ) coincides with that previously 
derived in where the WW approximation has not been considered. We see that the 
DVCS amplitude with longitudinally polarized virtual photon is also sensitive to the spin 
structure of the nucleon. In principle, this part of the amplitude can be extracted from the 
data through the angular, spin and Q dependence of the differential cross section . 
This may provide us with additional possibility to probe J q in the nucleon. 



In the case of axial transverse SPD F„ we obtain from Eq. 



the following result for 



the N = 1 moment: 



du u F, 



WW, 



u,0 



His}) 2 



-i 



du u H(u, + e 2 {^i(A 2 ) + ^ 2 (A 2 )} 



+ other structures 



(47) 



For simplicity we do not write all possible spinor structures which appear in the RHS ([47|). 
For the comparison with the sum rules of Ref. JlO| we need only coefficient in front of 
structure ((7^75)). We obtained that our result ( f4"TD is in agreement with Ref. flIP| . In the 
forward limit we reproduce from the Eq. (|47D the Efremov-Leader-Teryaev sum rule [25] for 
the polarized structure function grix). 

Discontinuities: Using general properties of the WW transformation ([|1|) and Eqs. ( |33| ) and 
( p5|) one obtains that F^ w and F^ w have discontinuities at the points x = ±£ . But using a 
certain symmetry of the equations (|33|) and (0), one can find that some combinations of the 
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distributions F ww and F ww are free of discontinuities. For example, using formulae ( flip 
one can see that the combination: 



WW, 



has no discontinuity at x = £. On the other hand, the 'dual' combination: 



F^ W (x,0+is±, P F p ww (x,0 



(4? 



(49) 



is free of the discontinuity at x = — £. As we shall see below, the cancellation of discontinuities 
in these particular combinations of the SPD's ensures the factorisation of the twist-3 DVCS 
amplitude on the nucleon. 



Pion pole and D-term in and F^ 

In the Wandzura-Wilczek approximation the new skewed distributions F^ and F^ (see Eq. (|7|) 
for definition) are fixed completely in terms of the twist- 2 SPD's H, E and H, E. These 
functions have some specific contributions which are invisible in the forward limit. In this 
section we consider how these special contributions are incorporated into the structure of 
the twist-3 distributions. 

Pion pole contribution: The skewed parton distribution E at small t is dominated by chiral 
contribution of the pion pole [EBL E3, [30, |31j of the form: 



Ag\M 2 1 (x\ „„ . 



(50) 



where gr^ is the axial charge of the nucleon and ^(w) is the pion distribution amplitude. In 
the WW relation ( |33"1) the twist- 2 SPD enters only in the combination: 



l + u— + e— 
<9« <9£ 



£(u,0 



(51) 



One can easily see that the contribution of the pion pole ( |50|) nullifies under action of the 
differential operator in Eq. fl5l|) . The only place where the pion pole contribution survives 



is the SPD F„, see Eq. 



In this way we come to the following simple results for the 



contribution of the pion pole to the twist-3 SPD's in the WW approximation: 



pWW,pion pole^ ^ 
p<WW,pion pole^ Q 



o, 

i3lM2 



2M 



1 (x 
-t + ml U, 



<0 



(52) 



We see an interesting result that the pion pole contribution to the twist-3 SPD's is expressed 
in terms of twist- 2 pion distributions amplitude. This observation is in nice agreement with 
the fact that in the WW approximation there exists no twist-3 pion distribution amplitude 
associated with vector- and axial-vector operators P2L . 
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D-term contribution: Now we discuss a specific contribution to the twist- 2 SPD's H and E 
which was found in Ref. [34|. This contribution completes the parametrisation of SPD's in 
terms of double distributions lEHl. This contribution has the form: 



x 



x 



E D - tc ^(x,0 = -d\j\ e(\x\<t) 



(53) 



where the function D(u) is odd in its argument. The specifics of the D-term contribution is 
that it is not 'generated' by double distributions. Numerically it can be rather large as it 
follows from estimates in the chiral quark-soliton model of the nucleon [Q. The importance 
of the D-term is especially emphasized by the fact that it contributes considerably to the 
DVCS amplitude but drops out in Ji's angular momentum sum rule @ because in the 
combination H + E D-term cancels. Now if we substitue the D-term contributions (1531) to 



our WW relations (|33|,|35|) we see that the result is: 



iWW D-term, 



w <0 



^WW D-term ^ = Q 



(54) 



Again we reproduced the result anticipated in Ref. |TTJ . 



DVCS amplitude 

Here we analyse the general structure of the DVCS amplitude on the nucleon to the twist-3 
accuracy. To this accuracy the amplitude receives contributions not only from the transverse 
polarization of the virtual photon (as it is at the twist- 2 level), but also from the longitudinal 
one. 

The amplitude with the transverse polarization of the virtual photon has the formQ (see 
Eq. ©) 



dx e-t(q) 



w 



F + (x,Z)C + (x,0 



-.v k\ 



P u A k ± 
W 



ieiF + (x,OC-(x,0 



(55) 



We see that in this case the amplitude to the 1/Q accuracy is fixed completely in terms of 
the same combinations of the twist- 2 SPD's. Such result for the DVCS amplitude to the 1/Q 
accuracy confirms the ad hoc prescription suggested in Ref. [0, [I8| . However the twist-3 
DVCS amplitude also receives contribution from the longitudinal polarization of the virtual 
photon (see Eq. 



**We neglected the last term in the expression (g) for the amplitude as it being contracted with polarization 
of the final real photon contributes only to the order 0(1/ Q 2 ) which is beyond our accuracy 
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e L (q)T^ cx 



P u A k 

3 ,±_ W 



X 



J^dx{F k (x,£)C + {x,S) -ie ±kp F p (x,Z)C-(x,Z)} . (56) 



This part of the amplitude depends on the 'transverse' parts of SPD's F^ ± and F fl± which 
are related to twist- 2 SPD's H, E and H, E only in Wandzura-Wilczek approximation. Since 
the coefficient functions C ± (x, £) have singularities at points x = ±£ the integral over x may 
diverge due to discontinuities of SPD's in WW approximation. However it is easy to see that 
near the point x = £ SPD's enter in the combination (081) and near the point x = — £ in the 



combination (f49|). These are exactly the combinations which are free of discontinuities at the 
corresponding point, therefore the amplitude (|56| ) is well defined in the WW approximation. 
The same cancellation of the discontinuities was observed in the case of DVCS on the pion 



target P"2"| , W . In Ref. [TEfl some general arguments why this phenomenon should persist to 



the case of DVCS on any target was given. Here we have checked by the explicit calculations 
this conjecture for the case of the nucleon target. 

Conclusions 

We have presented a simple and detailed derivation of the Wandzura-Wilczek relations for 
the twist-3 skewed parton distributions in the nucleon. These relations allow to express the 
DVCS amplitude on the nucleon to the twist-3 accuracy (1/Q) in terms of the twist-2 SPD's 
H, E and H, E. 

The Wandzura-Wilczek approximation corresponds to neglecting the nucleon matrix ele- 
ments of the twist-3 quark-gluon operators. It would be very interesting to check in various 
models of the nucleon structure whether the neglected matrix elements are indeed small. 
The theory of instanton vacuum predicts that in the case of polarized structure functions 
gr{x) and the corresponding matrix elements are parametrically small in the packing 
fraction of the instantons in the QCD vacuum, for details see |14], |I5| . Another interesting 
possibility to check the hypothesis behind the WW approximation is the model where SPD's 



are represented as an overlap of nucleon wave functions [37]. The quark-gluon correlations 
appear in this case as the overlap between the lowest Fock component of the nucleon wave 
function and that with an additional gluon. 

We have checked by an explicit calculation that although the SPD's in the WW approx- 
imation have discontinuities, the corresponding would-be divergences in the twist-3 DVCS 
amplitude are cancelled exactly and do not lead to problems with the factorization. Of 
course, this is not a proof of the factorization for the twist-3 DVCS amplitude. 

We rederive the new angular momentum sum rule of Ref. [|Kj] which allows us to relate 
twist-3 SPD's to the angular momentum fraction carried by quarks in the nucleon. Note that 
this sum rule is valid beyond the WW approximation. Since, in principle, the twist-3 SPD's 
can be extracted from the experimental data by considering Q, spin and angular dependencies 
of the DVCS observables, this opens a complementary approach to the measurement of the 
quark angular momentum distribution, compared to the original proposal by Ji H. 
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The twist-3 contributions are suppressed in the DVCS amplitude relative to the twist- 
2 one by one power of 1/Q only. Therefore, estimates of these corrections are extremely 
important for the possibility to extract twist- 2 SPD's from experimental data. The size 
of the twist-3 corrections increases with the momentum transfer squared t, because these 
corrections typically enter the amplitude as \f—tjQ. Their studies are hence, also helpful 
for the analysis of the t— dependence of twist-2 SPD. The latter task is mandatory for the 
very possibility to extract information about L q from the data because L q is related to the 
SPD's only at t — 0. This point lies outside the physical region of the DVCS process so that 
extrapolation is required. 
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APPENDIX: 



For the convenience of the reader who may wish to use our results for the calculation of 
DVCS observables to the twist-3 accuracy, we collect the resulting expressions for the DVCS 
amplitude to the accuracy 0(1/Q) in the Wandzura-Wilczek approximation. We present 
the results in the form which can be easily used for the analysis of DVCS observables. 



The DVCS amplitude on the nucleon to the order 0(1/Q) obtained in Refs. [^, [TT| has 
the form: 



dx 



:-<ru - 



nPF p {x,i)C+{x,i) 



(q + A^PY 1 



P u A k L 
~{Pq) 

(-g uk ) 



(Pi) 

ieiyF p {x,OC-{x,i) 



(Pq) 

where to the twist-3 accuracy: 

1 



P U A\ 
JPq) 



{F k (x, 0C + (x, - iei p F p (x, 0C~(x, £)} 



P 



ip + p') 



n 



A 



p — p 



-2£P + A_ 



o 2 



q = q 



A = — n 

4£ 



A 



(A.l) 



where p,p' are momenta of the initial and final nucleon and q, q' are momenta of the initial 
and final photon respectively. The light-like vectors n and n* are normalized as (n ■ n*) = 1. 
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Also we introduce the metric and totally antisymmetric tensors in the two dimensional 
transverse plane (£0123 = 



(-(T)± = ~f v + nV + n v n*\ = t^ aP n a P p 
The leading order coefficient functions are: 



(A.2) 



^,0 = 



± 



x — i + is x + £ — ie 



The skewed parton distributions and can be related to the twist-2 SPD's H, E, H and 
E with help of Wandzura-Wilczek relations: 

F? W (x,0 = ^((jj))E(x,0-^((l + ))(H + E)(x,0 + 



2£ 



+ 



J du G^u, £)W+(x, u, + ie ±ltk J du G k (u, £)W-(x, u, , 



1 // 75 



F™ w (x,0 = A,-((p)E(x,0-^((l + l 5 ))H(x,0 + 



2£ 



+ J du G^(u,£)W+(x,u,£) +ie ±fMk J ^du G k (u,£)W-(x,u,£) 
The following notations are used: 



((^}}(H + E)(u,0 + ^((^ 



2e 



«7 + » 



u _d_ + 
du d£. 



2£ M 

H + E)(u,0 



d , d 
U du- + ^ 



E(u,0~ 



1 + u d_ + ^d_ 

du dt; 



E(u,0 



AM 

•^((7 + 75)) 



du <9£ 



H(u,0 



rrr , + S 1 f „ / ^ 6 (u > x) „ . ^ 6 (u < x) ~] 

W ± (x,u,Z) = -|g(x>0 1 > -6{x<i) \_' ] 



6{u > x 



6{u < x) 



±-<0(x > -0— - 6>(x < -0— r 

21 v s; u + £ y s; M + f J 



The sandwiching between nucleon Dirac spinors is denoted by ((...)) = U(p') . . .U(p) and 
the quarks flavour dependence in the amplitude can be easily restored by the substitution: 

W)- E e}Ft(Ff). (A.3) 

f=u,d,s,... 
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